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ABSTRACT 


An  experimental  and  theoretical  investigation  of  bounded  turbulent 
mixing  with  chemical  reactions  is  described.  Configurations  which  are 
applicable  to  the  air-augmented  rocket  are  considered.  Experimental 
results  are  presented  for  a  mixing  apparatus  in  which  a  fuel-rich  O2-H2 
rocket  stream  was  mixed  with  a  secondary  air  stream  inside  a  conical 
duct.  An  integral  mixing  theory  is  described,  in  which  conservation 
equations  are  satisfied  across  the  duct.  The  inviscid  portions  of  the 
duct  flow  are  considered  one-dimensional,  and  the  mixing  zone  profiles 
are  assumed  to  exhibit  shape  similarity.  The  theory  is  extended  to  in¬ 
clude  the  downstream  regime  where  the  mixing  region  extends  across 
the  entire  duct.  Correlations  of  the  theory  with  incompressible  and 
reactive  compressible  mixing  experiments  indicate  that  the  theory  pro¬ 
vides  a  satisfactory  overall  representation  of  the  bounded  mixing 
process. 
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SECTION  I 
INTRODUCTION 


There  is  currently  considerable  interest  in  air -augmentation  as  a 
method  for  significantly  increasing  the  performance  of  fuel-rich  chemical 
rockets  during  atmospheric  flight.  Of  particular  interest  is  the  configura¬ 
tion  shown  in  Fig.  1,  which  consists  of  a  divergent  mixing  duct  connected 
to  an  annular  air  inlet  located  near  the  missile  base.  The  rocket  exhaust 
and  air  streams  mix  and  burn  inside  the  duct,  and  a  net  increase  in  thrust 
is  obtained  if  the  axial  pressure  force  on  the  mixing  duct  inner  surface  is 
larger  than  the  combined  inlet  and  external  drag.  The  weight  of  the  addi¬ 
tional  hardware  must  also  be  considered. 

Divergent  air -augmentation  configurations  were  first  investigated 
experimentally  in  1950  (Ref.  1)  for  the  case  where  the  secondary  stream 
stagnation  pressure  was  equal  to  the  ambient  pressure  at  the  duct  exit 
(corresponding  to  zero  flight  velocity).  More  recently,  divergent  air- 
augmentation  configurations  have  been  proposed  (Ref.  2}  for  use  at  high 
vehicle  velocities  with  substantial  ram  compression  of  the  secondary  air. 
The  performance  of  such  air-augmented  rockets  has  been  analyzed  by 
assuming  complete  mixing  and  burning  (or  some  fraction  thereof)  with 
arbitrary  pressure  distributions  along  the  mixing  duct  (Ref.  2),  Although 
this  type  of  one -dimensional  analysis  serves  to  point  out  the  overall  poten¬ 
tial  of  the  air -augmented  rocket,  it  is  incapable  of  predicting  the  duct  shape 
required  to  realize  an  assumed  pressure  distribution,  or  of  predicting  the 
performance  attainable  with  mixing  ducts  of  practical  shape.  A  perceptive 
theory  for  the  complex  process  of  bounded  mixing  and  burning  is  required 
for  realistic  evaluation  and  optimization  of  the  air-augmented  rocket. 

The  phenomenon  of  bounded  mixing,  with  or  without  burning,  is  also  en¬ 
countered  in  other  devices,  for  example,  the  supersonic  combustion  ram¬ 
jet  and  the  jet  pump.  For  a  jet  pump  with  a  constant  area  mixing  duct, 
one -dimensional  theory  satisfactorily  predicts  the  overall  performance; 
however,  the  length  of  duct  required  to  attain  the  performance  predicted 
by  one -dimensional  theory  can  be  predicted  only  by  a  detailed  analysis  of 
the  mixing  process. 

This  report  summarizes  the  results  of  an  experimental  and  theoreti¬ 
cal  investigation  of  the  bounded  coaxial  stream  mixing  problem.  The 
fundamental  objective  of  this  investigation  is  to  determine  (for  a  given 
central  stream,  duct  configuration,  and  secondary  stream  stagnation 
pressure)  the  induced  secondary  flow  and  duct  wall  pressure  distribution. 
Because  of  the  very  complex  nature  of  this  problem,  many  simplifications 
are  required  in  a  theoretical  treatment  of  the  overall  duct  flow.  Use  of 
integral  conservation  equations,  however,  provides  overall  results  which 
are  sufficiently  accurate  for  engineering  analysis,  even  though  the  detailed 
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mechanisms  are  over-simplified.  Emphasis  has  been  placed  on  relatively 
long  mixing  configurations  in  which  the  mixing  region  encompasses  most 
or  all  of  the  duct  cross  section  and  in  which  the  mixing  cannot  be  considered 
a  perturbation  on  the  inviscid  flow  fields  which  would  exist  without  mixing. 

The  bounded  mixing  problem  can  be  divided  into  three  distinct 
regimes  as  shown  in  Fig.  2.  In  the  first  regime,  turbulent  mixing 
occurs  between  the  inviscid  central  and  outer  streams.  In  the  second 
regime,  the  central  inviscid  stream  has  been  dissipated,  but  a  region  of 
inviscid  secondary  flow  exists  near  the  duct  wall.  The  third  regime 
occurs  when  the  mixing  zone  spreads  to  the  wall.  In  the  first  two  re¬ 
gimes,  the  mixing  is  essentially  free  turbulent  in  nature,  although  at 
variable  pressure,  and  free  turbulent  spreading  relations  may  be  applied. 
The  third  regime,  however,  requires  a  departure  from  the  techniques 
used  in  the  analysis  of  free  turbulent  mixing. 

The  mixing  zone  chemistry  is  assumed  to  be  either  frozen  or  in 
equilibrium;  that  is,  the  chemistry  is  considered  to  be  either  very  slow 
or  very  fast  with  respect  to  a  characteristic  time  for  the  flow. 


SECTION  II 

FORMULATION  OF  MIXING  THEORY 


Theoretical  treatments  of  free  turbulent  mixing  with  chemical  reac¬ 
tions  are  usually  limited  to  constant  pressure,  or  to  the  case  of  a  pre¬ 
scribed  axial  pressure  distribution  {e.  g. ,  Refs.  3,  4).  The  presence 
of  a  wall  in  proximity  to  the  mixing  region,  however,  causes  the  axial 
pressure  distribution  to  be  strongly  dependent  upon  the  mixing  process. 

The  bounded  mixing  problem  then  involves  simultaneous  solution  of  the 
inviscid  stream  flow  conditions  along  with  the  mixing  zone  conditions. 

The  integral  technique,  in  which  overall  conservation  equations  are  satis¬ 
fied  across  the  mixing  region,  is  used  in  this  analysis.  Use  of  the  integral 
method  is  advantageous  in  that  the  overall  results  (i.  e. ,  wall  pressure 
distribution)  are  relatively  insensitive  to  the  choice  of  mixing  zone  pro¬ 
files.  Similar  techniques  have  been  used  by  Mikhail  (Ref.  5)  and 
Yakovlevskiy  (Ref.  6)  for  incompressible  bounded  mixing. 


A  useful  relationship  between  the  velocity,  enthalpy,  and  composi¬ 
tion  distributions  can  be  obtained  by  considering  the  conservation  equa¬ 
tions  in  differential  form,  with  turbulent  Prandtl  and  Lewis  numbers  of 
unity  (Ref.  3): 


Momentum  equation: 

d  u 

pu  IT  +  pv 


d  u 
dr 


jr(ptT 


(Ju  \ 

jr) 


i  d  P 


(1) . 
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Energy  equation: 


pu 


d  H, 


dT^+  Pv  ^7  -  ■  -j-' 

Conservation  of  elemental  species: 


H , 


pu 


4=-  *  pv  4£-  -  r- 
ox  r  dr 


■  -t  (-« -if) 


(2) 


(3) 


Global  continuity  equation: 


■fr  (pu)  *  r  *  (pvr)  =  0 


(4) 


These  equations  are  the  laminar  boundary -layer  equations  with  the 
laminar  transport  properties  replaced  by  the  corresponding  effective 
turbulent  values.  Consider  the  case  of  constant  pressure  mixing  of  an 
initially  uniform  jet  with  a  surrounding  uniform  infinite  stream.  The 
dp/dx  term  in  Eq.  (1)  is  then  zero,  and  Eqs,  (1),  (2),  and  {3)  are  identi¬ 
cal  in  form.  A  linear  relation  is  obtained  between  the  variables  u,  Hn, 
and  C  : 


_  ~  Non  _  C  -  Ca  _  u  —  ua'  -• 

H0J  —  Hoa  Cj  —  Ca  Uj  ~  ua  (^) 


The  stagnation  enthalpy,  H0.  is  defined  to  include  the  chemical  heats 
of  formation.  For  constant  pressure  mixing,  H0  and  C  are  related  to 
velocity  by  Eq.  (5),  and  solution  of  momentum  and  global  continuity  equa¬ 
tions  completely  defines  the  flow  field.  For  application  to  an  integral 
method,  Eq.  (5)  is  also  assumed  to  be  valid  for  mixing  with  moderate 
axial  pressure  gradients.  Note,  however,  that  the  inviscid  stream  refer¬ 
ence  velocities  are  now  pressure  dependent. 

The  additional  assumptions  made  in  the  integral  treatment  for  bounded 
mixing  are  as  follows: 

1.  -  The  inviscid  streams  are  one-dimensional  and  isentropic. 

2.  The  static  pressure  is  constant  across  any  duct  section. 

3.  Wall  viscous  effects  on  the  velocity  profiles  are  neglected.  c 
However,  for  very  long  mixing  ducts,  the  wall  boundary- 
layer  effect  may  be  included  as  a  wall  displacement. 

4.  The  velocity  profile  at  the  initiation  of  mixing  is  a  step 
function. 

5.  If  the  central  stream  is  supersonic  and  not  correctly  expanded, 
a  one-dimensional  isentropic  adjustment  process  is  assumed  in 
which  the  central  and  outer  streams  expand  or  contract  until 
pressure  equilibrium  is  reached  and  the  combined  stream  areas 
equal  the  initial  duct  area. 


\ 
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2.1  FIRST  REGIME 


The  integral  conservation  equations  in  axisymmetric  coordinates 
{Fig.  3)  are: 


Axial  momentum  across  the  duct: 
d 


f  pu2  rdr  =  -  —L-  f  rdr 


Steady  flow  continuity: 

/  Purdr  =  0 

o 

These  equations  become  for  the  three  regions  of  the  flow  (central 
inviscid  stream,  mixing  zone,  and  outer  inviscid  stream): 


(6) 


(7) 


3 

dx 


{'W  4- +  '■■■■■  *"y]  *  /  -  -  -Hr  Jf 

4-+  <>•»*  +  h 'I  *  /+l  o.*}  -  o 


(6a) 

(7a) 


The  inviscid  stream  densities  and  velocities  (Pa,  Pj,  ua,  uj)  are  isentropic 
functions  only  of  p. 


2.1.1  Mixing  Zone  Profiles 


For  free  mixing  with  moderate  axial  pressure  gradients,  it  is  well 
established  that  experimental  mixing  zone  velocity  profiles  exhibit  shape 
similarity  at  various  axial  stations: 


U  -.11 


a 


u 


J 


Various  profile  shape  equations  are  available  which  correlate  reason 
ably  well  with  experimental  velocity  distributions.  The  integral  method 
is  not  sensitive  to  the  profile  chosen  so  long  as  it  is  a  reasonable  repre¬ 
sentation  of  the  experimental  shape.  A  cosine  profile  was  used  because 
this  profile  reaches  the  inviscid  stream  velocities  at  finite  radial 
distances, 

U  —  11,  _  l 

Uj  -  u.  2 


[i +  cos *  (-v^)] 

For  frozen  mixing,  Eq.  (5)  may  be  interpreted  as 

hr,  —  hr,.  C  —  C.  U  —  U_ 


h  o j  “  h0  H 


C  -c. 
Cj-Cfl 


(8) 


(5a) 
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where  h0  is  the  sensible  enthalpy.  Because  there  is  no  chemical  reac¬ 
tion,  C  may  be  taken  as  the  mass  fraction  of  any  specie.  The  mixing 
zone  specific  heat  and  gas  constant  are  evaluated  as  mass  averages  of 
the  constituent  values. 


For  mixing  with  chemical  reactions,  Eq.  (5)  is  defined  as 


*0  8  _ 


H0J  —  H0 


=  -  C 

Cj  -  C. 


(5b) 


where  C  is  the  mass  fraction  of  elements  originating  in  the  central 
stream  (thus  Ca  =  0  and  Cj  =  1).  This  composition  parameter  is  con¬ 
venient  when  the  inviscid  streams  consist  of  several  elemental  species. 
The  assumption  of  equilibrium  mixing  zone  chemistry,  along  with  the 
energy  equation  and  the  perfect  gas  law,  then  gives  the  mixing  zone 
temperature  and  density  for  any  velocity.  This  treatment  of  the  mixing 
zone  implies  that,  at  any  point  on  the  non-dimensional  profile,  the  mass 
fractions  of  the  elements  are  the  same  as  for  frozen  mixing  and  that 
chemical  reaction  shifts  only  the  local  temperature  and  gas  properties 
from  the  frozen  values. 


To  simplify  the  computation  procedure  for  the  equilibrium  chemistry 
case,  temperature  and  gas  property  functions  of  C  were  computed  at 
representative  values  of  p,  uj,  and  ua,  using  the  method  presented  in 
Appendix  I.  The  distribution  of  total  temperature  was  then  computed, 
and  the  gas  property  and  total  temperature  functions  of  C  were  taken  as 
constant  throughout  a  computation.  By  using  this  technique,  the  chem¬ 
istry  enters  the  mixing  problem  only  through  the  gas  property  and  total 
temperature  functions  of  C  ;  a  typical  set  of  these  functions  is  shown  in 
Figs.  4,  5,  and  6. 

2.1.2  Turbulent  Spreading  Rate 


The  mixing  zone  profiles  of  p  and  u  are  now  known  with  respect  to 
dimensionless  position  in  the  mixing  zone.  A  relation  expressing  the 
rate  of  growth,  db/dx  ,  of  the  mixing  zone,  is  required  to  allow  solution 
of  Eqs.  (6a)  and  (7a)  for  p(x)  and  r,(x).  Abramovich  (Ref.  7)  proposed 
the  following  relationship  for  free  turbulent  mixing  with  large  density 
gradients: 


■lb 

ilx 


(*W  “  Pa°»J)  J  Pd  (-T-'  ) 


u 


pu  d  j 


f  r  “ 


*■>]' 


(9) 


where  kj  is  one -half  the  incompressible  mixing  zone  spreading  rate  for 
ua  =  0.  The  value  for  kj  is  established  as  about  0.  13  (Ref.  8);  however, 
the  experimental  determination  of  k|  from  mixing  zone  velocity  measure¬ 
ments  is  very  difficult  because  the  outer  edges  of  the  mixing  zone  are 
not  easily  located.  It  is  also  probable  that  the  analytical  velocity  pro¬ 
file  used  in  a  particular  theory  is  not  an  exact  representation  near  the 


5 


AEDC-TR-65-4 


low  velocity  boundary;  therefore,  a  better  procedure  for  determining 
L[  is  to  correlate  the  particular  theory  with  incompressible  experiments 
to  give  the  best  overall  representation  of  the  flow. 

2.1.3  Elemental  Species  Conservation 

The  relation  between  mixing  zone  velocity  and  concentration  (Eq.  5) 
was  deduced  from  the  conservation  equations  for  the  case  of  constant 
pressure  mixing.  The  same  relation  was  assumed  to  be  valid  for  vari¬ 
able  pressure  mixing,  and  the  validity  of  this  assumption  is  checked 
throughout  the  mixing  theory  computations.  Continuity  of  the  central 
stream  elemental  species  is  checked  by  evaluating  the  function, 

Q  =  ~ ~  f  puC  rdr  (lt 

If  Q  does  not  deviate  from  unity,  then  elemental  species  continuity  is 
satisfied,  and  use  of  Eq.  (5)  is  justified,  at  least  for  use  in  an  integral 
method. 


2.2  SECOND  REGIME 


The  second  regime  is  treated  in  the  same  manner  as  the  first, 
except  that  the  mixing  zone  orientation,  rj,  is  replaced  as  a  variable 
by  the  centerline  velocity,  un.  Equations  (6)  and  (7)  become 


The  profile  shape  equation  becomes 


The  turbulent  spreading  rate  equation  becomes 
^  _  kll  (Pc  “c*  -  P.  “a3)  „/1  P< I  (l~) 

dx  "  2  [jVM-ir)]7 


(8a) 


(9  a) 


The  value  of  the  second  regime  spreading  parameter  is  reported  as  0.  11 
(Ref.  8),  but  again,  it  should  be  evaluated  by  correlating  the  theory  with 
incompressible  experiments. 
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2.3  THIRD  REGIME 


Mikhail  (Ref.  5)  and  Yakovlevskiy  (Ref.  6)  have  shown  that,  for  in¬ 
compressible  flow,  the  velocity  profiles  in  the  third  regime  have  essen-  . 
tially  the  same  shape  as  free  mixing  profiles.  Of  course,  very  far 
downstream  the  wall  effects  predominate,  and  the  profile  loses  the  free 
mixing  shape;  however,  for  a  few  diameters  downstream,  the  wall  shear 
is  small  compared  to  the  free  turbulent  shear,  and  the  assumption  of 
profile  shape  similarity  is  justified.  The  concept  of  mixing  zone  growth 
rate  is  no  longer  valid  in  the  third  regime;  therefore,  the  turbulent 
spreading  relation  is  replaced  by  an  additional  momentum  equation  written 
from  the  centerline  to  the  duct  half  radius  (Fig.  7): 


d_ 

d  x 


/ 


0 


P  «*’  rdr  ~  UM 


_d_ 
d x 


o 


(11) 


where  the  subscript  M  refers  to  the  duct  half  radius  position.  For  a  con¬ 
siderable  distance  downstream,  the  turbulent  shear  at  the  profile  mid¬ 
point  will  be  free  turbulent  in  nature,  and  may  be  expressed  as 

'  =  W  “If |m  (12) 

where  (  is  the  free  turbulent  eddy  viscosity.  The  Prandti  formulation 
for  e  is 


f  —  k  b  (llmax 


Um  in  * 


(13) 


where  e  is  assumed  constant  across  the  mixing  zone.  It  is  well  known  that 
t  is  not  constant  laterally  for  incompressible  flow  and  is  certainly  not  for 
mixing  with  large  density  gradients.  The  procedure  used  in  this  study 
was  to  assume  that  Eq.  (13)  is  functionally  valid  along  the  velocity  pro¬ 
file  midpoint  and  may  be  used  to  extrapolate  the  half  radius  shear  stress, 

7m  ,  which  is  evaluated  by  solving  Eq.  (11)  at  the  end  of  the  second  regime. 
The  half- radius  shear  stress  is  then  at  any  station,  replacing  b  by  t„, 

rM  *  Pm  _  Uw)1 


and  the  relation  used  to  relate 

rM  =  FM 


to  7m  is 

l°M  (“c  ~  O* 
PM  (ue  -  “w)J 


(14) 


where  the  barred  variables  are  those  evaluated  at  the  end  of  the  second 
regime. 


The  three  conservation  equations  (6),  (7),  and  (11)  may  now  be  solved 
for  the  unknowns,  uc,  uw,  and  p.  The  integral  conservation  equations 
have  been  transformed  to  a  form  suitable  for  numerical  computations 
(Ref.  9)  and  are  solved  with  a  high-speed  digital  computer. 
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SECTION  Ml 

EXPERIMENTAL  PROGRAM 


The  RTF  Propulsion  Research  Cell  (R-1B)  (Fig.  8)  was  designed  to 
investigate  the  bounded  mixing  and  burning  process.  A  small  scale  rocket 
was  used  to  provide  the  primary  stream,  and  atmospheric  air  was  used 
for  the  secondary  stream.  A  conical  mixing  duct  was  used.  The  cell 
was  connected  to  the  RTF  exhauster  system,  which  maintained  the  back 
pressure  on  the  mixing  duct  at  approximately  2.  5  psia.  The  principal 
dimensions  and  operating  parameters  for  the  apparatus  are  given  in 
Table  I. 

In  addition  to  the  usual  measurements  of  rocket  and  air  stream 
parameters,  the  following  parameters  were  measured; 

1.  Mixing  duct  static  pressure  distribution, 

2.  Pitot  pressure  distribution  across  mixing  duct 
exit  plane,  and 

3.  Gas  composition  distribution  across  mixing  duct 
exit  plane. 


3.1  ROCKET  ENGINE 

A  water-cooled,  oxygen -hydrogen  rocket  engine  was  used  to  gener¬ 
ate  a  high  temperature  stream  of  fuel-rich  exhaust  gases.  The  character¬ 
istic  length,  L*  (chamber  volume/throat  area),  of  the  rocket  chamber  was 
approximately  35  in. 

The  oxygen-hydrogen  propellant  combination  was  chosen  because  of 
the  relatively  simple  exhaust  gas  composition  which  lends  itself  readily 
to  gas  analysis.  Gaseous  propellants  were  chosen  to  achieve  a  high  com¬ 
bustion  efficiency;  liquid-propellant  rocket  engines  of  this  scale  normally 
operate  at  much  lower  combustion  efficiency  than  a  comparable  large 
booster  engine.  A  rocket  engine  having  high  combustion  efficiency  is 
necessary  for  accurate  experimental  evaluation  of  rocket  thrust  aug¬ 
mentation  configurations  in  which  afterburning  of  the  rocket  exhaust 
occurs.  If  the  rocket  is  inefficient,  an  unrealistically  large  amount  of 
unburned  fuel  is  available  for  afterburning.  This  excessive  afterburning 
can  lead  to  erroneous  conclusions  about  the  performance  of  the  augmented 
rocket  as  compared  to  that  of  the  basic  rocket. 

The  theoretical  combustion  temperature  of  the  gaseous  02-H2  engine 
is  approximately  5300°R.  The  characteristic  velocity,  c*  (defined  as 
chamber  pressure  times  throat  area  divided  by  propellant  flow),  is 
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normally  used  as  a  measure  of  combustion  efficiency.  For  the  25  fir¬ 
ings  discussed  in  this  report,  the  experimental  value  of  c*  averaged 
97  percent  of  theoretical  c*.  Most  of  the  deviation  from  ideal  c*  results 
from  heat  transfer  to  the  cooling  water. 


3.2  INSTRUMENTATION 


3.2.1  Wall  Static  Pressure 


Wall  static  pressures  throughout  the  apparatus  were  measured  with 
mercury  manometers,  referenced  to  atmosphere  and  recorded  photo¬ 
graphically.  Mixing  duct  static  pressures  were  measured  with  two  rows 
of  taps  90  deg  apart;  the  axial  tap  spacing  was  2  in. 

A  vacuum  check  was  made  prior  to  and  following  each  test  to  detect 
any  leaks  in  the  manometer  system.  The  scales  on  the  manometers 
were  subdivided  into  0.  1-in.  increments,  and  the  resulting  pressures 
were  read  with  an  estimated  precision  of  ±0.  05  psi. 

3.2.2  Air  Flow  Rata 

The  secondary  airflow  was  measured  with  a  choked  venturi,  having 
a  circular  arc  inlet  contour  and  a  conical  diffusion  section.  For  the  range 
of  throat  Reynolds  number  encountered,  the  nozzle  flow  coefficient  was 
greater  than  0.  99  (Ref.  10);  therefore  the  coefficient  was  assumed  to  be 
unity.  The  stagnation  pressure  of  the  venturi  flow  was  measured  with 
mercury  manometers,  and  the  stagnation  temperature  was  measured  with 
an  immersion-type  thermocouple. 

3.2.3  Propellent  Flow  Rates 

The  propellant  flows  were  also  measured  with  choked  venturis.  Pro¬ 
pellant  temperatures  were  measured  with  immersion-type  thermocouples 
which  were  located  downstream  from  the  venturis  to  avoid  distorting  the 
inlet  flow.  The  nozzle  flow  coefficients  were  again,  for  the  range  of 
throat  Reynolds  number  encountered,  greater  than  O’.  99  (Ref.  10). 

3.2.4  Rocket  System  Pressures 

Rocket  chamber  pressure  and  the  propellant  metering  nozzle  inlet 
pressures  were  measured  with  strain-gage -type  transducers  and  recorded 
on  direct -inking,  null-balance  potentiometers.  The  systems  were  periodic¬ 
ally  calibrated  against  a  secondary  standard  to  check  for  non-linearity  and 
for  absolute  level.  The  calibrations  were  checked  prior  to  each  test  period 
by  applying  a  fixed  resistance  to  the  system  to  obtain  a  full-scale  deflection. 
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During  shakedown  of  the  apparatus,  the  propellant  pressures  downstream 
from  the  metering  nozzles  were  measured  with  similar  instrumentation. 
These  pressure  measurements  were  discontinued  after  it  was  found  that 
the  pressure  drops  across  the  venturis  were  sufficient  to  ensure  choked 
flow  at  the  nozzle  throats. 

3.2.5  Temperatures 

Air  stream  temperature  was  measured  with  an  immersion-type 
thermocouple  and  recorded  on  a  multi-point,  null-balance  potentiometer. 

The  propellant  temperatures  were  measured  with  immersion-type 
thermocouples  and  recorded  on  a  light-beam  oscillograph.  The  thermo¬ 
couple  system  was  calibrated  by  applying  known  voltage  from  a  standard 
cell  and  recording  the  galvanometer  deflection.  The  thermocouples  were 
referenced  to  control  room  temperature  which  was  measured  with  a 
mercury  bulb  thermometer. 

3.2.6  Survey  Rake 

The  mixing  duct  exit  plane  surveys  were  made  with  a  water-cooled, 
thirteen- probe  rake  having  seven  pitot  pressure  probes  and  six  gas  sam-' 
pling  probes  (Fig.  9).  The  rake  was  installed  with  the  probe  tips  about 
1/4  in.  downstream  from  the  mixing  duct  exit;  a  complete  survey  of  the 
exit  plane  flow  was  accomplished  by  installing  the  rake  during  different 
tests  with  various  spacers  at  the  rake  mounting  flange.  Details  of  the 
probe  tip  construction  are  shown  in  Fig.  9b. 

The  gas  sampling  probes  were  similar  to  those  used  by  Rhodes,  et  al 
(Ref.  11)  in  their  investigations  of  shock-induced  combustion  of  hydrogen- 
air  mixtures.  A  sapphire  watch  bearing  (fused  aluminum -oxide)  was  im¬ 
bedded  in  the  probe  tip  to  provide  an  orifice  capable  of  withstanding  the 
high  stream  temperatures.  The  probe  was  designed  so  that  the  expansion 
of  the  sample  flow  inside  the  probe,  and  the  subsequent  cooling  effect  of 
the  inner  probe  walls,  would  cause  quenching  of  the  major  chemical  reac¬ 
tions.  The  gas  samples  were  taken  in  evacuated  bottles  and  were  later 
analyzed  for  N2*  O 2 ,  and  H2  with  a  gas  chromatograph. 

The  rake  pitot  pressures  were  measured  with  mercury  manometers 
referenced  to  atmosphere.  Near  the  duct  centerline,  however,  the  pres¬ 
sures  were  too  high  to  be  measured  with  the  available  100 -in.  manometers, 
and  dial -type  bourdon-tube  gages  were  used.  These  gages  were  graduated 
in  0.  2-psi  increments  and  could  be  read  with  an  estimated  precision  of 
±0.  1  psi. 
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3.3  EXPERIMENTAL  PROCEDURE 

'■  The  desired  airflow  was  set  by  adjusting  a  throttle  valve  upstream 
of  the  venturi.  The  rocket  was  then  fired  for  a  nominal  duration  of 
30  sec.  Photographs  of  the  manometer  boards  were  taken  after  stabil¬ 
ization.  Gas  samples  were  taken  during  a  20 -sec  period  following 
rocket  stabilization. 


SECTION  IV 

RESULTS  AND  DISCUSSION 


4.1  INCOMPRESSIBLE  BOUNDED  MIXING 

The  mixing  theory  has  been  correlated  with  the  low  velocity  jet 
pump  experiment  of  Mikhail  (Ref.  5)  to  evaluate  the  theory  for  incom¬ 
pressible  flow  and  to  determine  the  incompressible  values  for  the  tur¬ 
bulent  spreading  parameters.  Mikhail's  apparatus  is  shown  in  Fig.  10, 
along  with  the  principal  experimental  parameters.  The  experimental 
axial  duct  pressure  distribution  is  shown  in  Fig.  11.  Theoretical  pres¬ 
sure  distributions  are  also  shown  for  two  cases:  (1)  wall  boundary-layer 
effect  neglected  and  (2)  wall  boundary -layer  effect  considered  as  an  esti¬ 
mated  linear  displacement  thickness.  Both  theoretical  curves  were  com¬ 
puted  with  kj  =  k[i  =  o.ll  (compared  to  0.  13  and  0.  11  reported  in  Ref.  8 
for  free  mixing).  The  theory  predicts  the  inflected  pressure  distribution, 
which  is  characteristic  for  the  low  velocity  jet  pump.  Theoretical  and 
experimental  axial  distributions  of  centerline  velocity  are  shown  in 
Fig.  12.  The  theoretical  distribution  agrees  well  with  the  experiment  in 
the  first  and  second  mixing  regimes,  and  the  agreement  is  satisfactory 
for  about  three  duct  diameters  downstream  into  the  third  regime.  Further 
downstream,  the  theory  does  not  predict  as  much  centerline  velocity  decay 
as  was  observed  experimentally.  This  indicates  that  the  flow  is  no  longer 
free  turbulent  in  nature  but  that  the  wall  effects  are  beginning  to  predomi¬ 
nate. 


4.2  02 -H2  ROCKET  EXPERIMENTS 

For  the  experiments  which  have  been  conducted  during  this  investiga¬ 
tion,  the  non-one-dimensional  effects  in  the  jet  plume  from  the  conical 
rocket  nozzle  cannot  be  neglected  when  applying  the  mixing  theory.  The 
rocket  exhaust  jet  expands  from  the  nozzle  and  causes  the  outer  air  stream 
to  choke  at  some  distance,  x*  ,  downstream  from  the  initial  section 
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(Fig.  13a).  The  outer  stream  is  supersonic  downstream  from  the  chok¬ 
ing  section.  Because  the  plume  is  not  one -dimensional,  the  jet  cross- 
sectional  area  is  larger  than  if  the  expansion  were  actually  one- 
dimensional.  The  secondary  flow  is  directly  proportional  to  the  secondary 
stream  choking  area,  Aa*  ,  and  the  following  procedure  was  used  to  deter¬ 
mine  that  area.  First,  the  streams  are  assumed  to  be  inviscid  up  to  the 
choking  station,  and  the  plume  shape  is  assumed  to  be  adequately  repre¬ 
sented  by  a  plume  issuing  from  the  rocket  into  an  ambient  medium  at  p*  . 
The  resulting  plume  area  distribution,  as  computed  by  the  method  of 
characteristics,  along  with  the  duct  area  distribution,  gives  a  distribu¬ 
tion  of  outer  stream  area.  The  minimum  area  is  taken  to  be  Aa  ,  and  the 
area  distribution  up  to  the  choke  point  defines  the  outer  stream  pressure 
distribution. 

The  mixing  theory  is  patched  onto  the  inviscid  plume  analysis  at  the 
choke  point,  and  the  difference  between  the  plume  cross-sectional  area 
and  one-dimensional  area  is  treated  as  a  small  parallel  wall  shift 
(Fig.  13b).  This  approximate  treatment  of  the  non-uniform  central  jet 
flow  must  be  regarded  as  tentative  and  is  inadequate  for  a  treatment  of 
the  flow  where  L  is  approximately  equal  to  x*.  For  long  mixing  ducts, 
where  L  is  much  larger  than  x*,  the  downstream  flow  is  relatively 
insensitive  to  the  non-uniform  initial  conditions,  and  the  approximate 
treatment  is  sufficiently  accurate  for  engineering  analysis. 

The  induced  air-rocket  mass  ratio,  wa/wj ,  for  the  O2-H2  rocket 
configuration  is  shown  in  Fig.  14,  plotted  against  the  stagnation  pressure 
ratio,  poa/poj  ■  The  theoretical  curve,  calculated  with  the  inviscid  plume 
analysis,  is  generally  higher  than  the  experimental  points,  with  less  than 
5  percent  deviation  over  most  of  the  experimental  range  of  poa/Poj  ■ 

Typical  experimental  mixing  duct  static  pressure  distributions  are 
shown  in  Figs.  15a  and  b  for  wa/wj  of  5.  3  and  3.  6,  respectively.  The 
theoretical  pressure  distributions  were  calculated  up  to  x*  by  the  inviscid 
plume  analysis,  and  downstream  of  x*  by  the  mixing  theory  with  the 
assumption  of  equilibrium  mixing  zone  chemistry  and  k[  =  0.10.  The 
theory  predicts  that  a  small  core  of  the  inviscid  rocket  stream  exists  at 
the  duct  exit;  therefore,  only  the  first  mixing  regime  is  encountered  with 
this  configuration.  The  experimental  pressure  distribution  reflects  the 
tendency  of  the  jet  to  alternately  expand  and  contract  because  of  the  non- 
one-dimensional  effects;  the  theory  represents  an  averaged  downstream 
distribution. 

The  results  of  a  pitot  pressure  survey  across  the  mixing  duct  exit 
plane  are  shown  in  Figs.  16a  and  b  for  the  same  experimental  conditions 
as  for  the  wall  pressure  distributions  in  Figs.  15a  and  b.  The  experi¬ 
mental  points  represent  the  results  of  several  tests  at  each  operating 
condition.  The  experimental  distributions  were  plotted  radially  from  the 
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centerline  of  the  profiles,  which  was  displaced  about  1/4  in.  from  the 
duct  centerline.  The  symmetry  of  the  experimental  profiles  can  be  seen 
by  noting  that  the  flagged  symbols,  which  represent  measurements  on  the 
opposite  side  of  the  duct,  generally  fall  along  the  same  distribution  as  the 
unflagged  symbols.  The  theoretical  distribution  in  Fig.  16a  is  seen  to 
agree  satisfactorily  with  the  experimental  distribution;  however,  the 
experimental  points  near  the  profile  centerline  in  Fig.  16b  fall  consider¬ 
ably  below  the  theory.  It  is  probable  that  this  deviation  results  from  sub¬ 
stantially  larger  stagnation  pressure  losses  in  the  inviscid  rocket  stream 
for  the  low  secondary  pressure  case  as  compared  to  the  higher  pressure 
case.  For  the  lower  pressure,  the  rocket  stream  initially  plumes  to  a 
greater  extent,  and  consequently,  the  downstream  shock  system  in  the 
inviscid  jet  is  stronger.  It  is  assumed  in  the  theory  that  these  inviscid 
stream  losses  are  negligible. 

The  results  of  exit  plane  gas  composition  surveys  are  shown  in 
Figs.  17a  and  b  for  the  same  experimental  conditions  as  for  the  data  shown 
in  Figs.  15  and  16.  'The  radial  distributions  of  N2.  O2,  and  H2  agree 
reasonably  well  with  those  predicted  by  the  theory,  with  the  assumption  of 
equilibrium  mixing  zone  chemistry.  The  experimental  points  shown  repre¬ 
sent  the  results  of  several  tests  at  each  operating  condition.  As  with  the 
pitot  pressure  data,  the  flagged  symbols  represent  samples  taken  on  the 
side  of  the  profile  centerline  opposite  to  that  for  the  unflagged  symbols. 

The  thrust  performance  of  this  mixing  duct  configuration  is  shown  in 
Fig.  18,  plotted  as  the  duct  thrust  integral,  Fd  =  /  pdA  ,  divided  by  the 
theoretical  rocket  vacuum  thrust,  Fj .  The  experimental  values  of  Fn  were 
obtained  by  graphical  integration  of  the  duct  pressure  distributions.  The 
duct  thrust  predicted  by  the  mixing  theory  with  the  initial  plume  analysis  is 
seen  to  agree  with  the  experimental  results  to  within  a  few  percent  over 
the  experimental  range  of  secondary  stream  stagnation  pressure.  The 
flagged  symbols  in  Fig.  18,  which  fall  consistently  below  the  other  points, 
represent  thrust  results  obtained  during  one  testing  period,  and  their 
validity  is  subject  to  question.  Careful  examination  of  all  the  experimental 
parameters  measured  during  that  testing  period,  however,  failed  to  yield 
a  definite  reason  for  the  experimental  deviation.  It  is  possible  that  the 
reference  pressure  for  the  manometers  on  which  the  duct  static  pressures 
were  measured  was  slightly  different  from  atmospheric  because  of  an 
obstruction  in  a  line. 

The  central  stream  species  parameter,  Q  ,  did  not  vary  more  than 
one  percent  from  unity  in  any  of  the  theoretical  computations  for  this 
configuration. 

Various  values  of  ki  were  used  to  correlate  the  theory  with  these 
experiments.  The  best  overall  correlation  was  obtained  with  kj  =  0.10, 
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as  compared  with  k|  =  0.11  used  to  correlate  the  incompressible  experi¬ 
ments.  If  the  Abramovich  turbulent  spreading  relation  were  exact,  the 
value  of  kl  should  also  have  been  0.  11  for  the  O2-H2  experiments.  The 
Abramovich  relation  is  thus  not  exact  but  predicts  a  slightly  excessive 
rate  of  mixing  for  the  case  of  very  large  density  gradients  in  the  mixing 
zone. 

The  mixing  zone  chemistry  appears  to  be  nearly  in  equilibrium  for 
these  experimental  conditions,  based  on  the  overall  correlation  between 
theory  and  experiment.  The  mixing  duct  static  pressure  distributions, 
and  radial  exit  plane  distributions  of  composition  and  pitot  pressure, 
are  satisfactorily  predicted  with  the  assumption  of  equilibrium  chemis¬ 
try  incorporated  in  the  mixing  theory. 


4.3  EFFECT  OF  PRESSURE  GRADIENTS 


The  turbulent  spreading  rate,  db/dx  ,  calculated  from  Eq.  9  does  not 
account  for  the  effects  of  axial  pressure  gradients.  It  is  well  known  that 
the  rate  of  growth  of  a  turbulent  wall  boundary  layer  is  appreciably 
affected  by  axial  pressure  gradients,  with  negative  pressure  gradients 
tending  to  decrease  the  rate  of  growth.  By  analogy  to  the  wall  boundary 
layer,  an  expression  was  derived  which  shows  that  even  small  axial 
pressure  gradients  can  significantly  influence  the  turbulent  spreading 
rate  when  the  width  of  the  mixing  zone  is  large.  By  writing  the  momen¬ 
tum  integral  equation  for  a  wall  boundary  layer,  with  the  inviscid  stream 
denoted  by  the  subscript  «  , 

A  6  Cfw  H  ■*  2  a  d  p 

dx  2  />»  u„  _  d* 

where  6  is  the  momentum  thickness  and  H  the  shape  parameter.  If  this 
equation  is  applied  to  the  free  turbulent  layer,  with  the  central  inviscid 
stream  as  the  reference  stream,  then 


d  0  _  _cj_  /8_ \  H  +  2  b  rip 
dx  2  \  h  /  pj  uj2  dx 

The  free  turbulent  friction  coefficient,  cf ,  is  a  function  of  the  local  flow 
conditions  and  not  of  dp/dx;  therefore,  the  rate  of  growth  of  the  turbulent 
zone  can  be  significantly  different  from  the  constant  pressure  rate  only 


if  the  quantity 


(4) 


H  +  2 
Pi  ui2 


dp 

dx 


i 


is  significant  with  respect  to  cf/2  , 


which  can  occur  if  either  b  or  dp/dx  is  large.  If  the  constant  pressure 
spreading  rate  is  denoted  by  the  subscript  o  ,  then 


db 

db 

re*  dpi 

L  2  \  b  /  Pj  Uj  rix  J 

rix 

dx 

o  -LL 

2 
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or 


db 


1  [1 +  2  {{ 

\  (H  +  2)  L  dp 

L  L  ^ h 

/  Cf  pj  UjJ  dx 

(15) 


If  Eq.  (15)  is  applied  to  the  high  velocity  half  of  the  mixing  zone,  the 
midpoint  friction  coefficient,  cfM,  is  used.  The  shape  parameter,  H  , 
and  momentum  thickness,  6 ,  are  evaluated  for  the  high  velocity  half  of 
the  layer.  These  quantitues  were  computed  for  a  mixing  region  typical 


for  the  02 -H2  rocket  configuration,  and  the  quantity 


(4)  -ii*-] 

V  b  /  CfM 


was 


found  to  be  about  10.  For  the  pressure  gradients  encountered  in  the 
02 -H2  rocket  experiments,  the  value  of  db/dx  was  only  slightly  different 
from  db/dx|  o  ,  and  the  pressure  gradient  effect  was  neglected.  The  pres¬ 
sure  gradient  effect  was  also  neglected  in  the  incompressible  jet  pump 
correlations.  This  tentative  analysis  does  point  out,  however,  that 
Eq.  (9)  can  be  considerably  in  error  for  mixing  configurations  which  have 
large  axial  pressure  gradients  or  large  mixing  zone  widths.  A  more  rigor¬ 
ous  treatment  of  the  pressure  gradient  effect  on  mixing  rates  would  be 
required  to  treat  such  configurations. 


SECTION  V 
CONCLUSIONS 


The  following  conclusions  can  be  drawn  from  the  results  of  this  investi¬ 
gation  of  bounded  turbulent  mixing  with  chemical  reaction: 

1.  The  integral  mixing  theory  provides  a  satisfactory  overall 
representation  of  the  bounded  turbulent  mixing  process. 

2.  The  assumption  that  the  elemental  species  concentration 
is  linearly  related  to  the  mixing  zone  velocity  appears  to 
intrinsically  satisfy  the  requirement  of  elemental  species 
conservation. 

3.  The  Abramovich  turbulent  spreading  relation  gives  results 
which  are  approximately  correct.  The  mixing  rates  for  a 
very  low  density  rocket  stream  are  satisfactorily  predicted, 
considering  the  current  state  of  knowledge  on  turbulent 
mixing. 

4.  Correlation  of  theoretical  and  experimental  results  for  the 
O2-  H2  rocket  configuration  indicate  that  the  mixing  zone 
reactions  are  nearly  in  equilibrium  at  the  mixing  duct  exit. 
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Further  investigation  is  required  on  the  following  problems  associated 
with  bounded  turbulent  mixing:  (1)  the  effect  of  pressure  gradients  on  tur¬ 
bulent  spreading  rates,  and  (2)  the  effect  of  non -uniformities  in  the  inviscid 
streams  on  the  overall  duct  flow. 
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APPENDIX  I 

EQUILIBRIUM  MIXING  ZONE  CHEMISTRY 


The  mixing  zone  composition  (elements)  varies  between  the  com¬ 
position  of  the  unmixed  central  stream  and  the  composition  of  the  outer 
air  stream;  thus,  the  elemental  fuel-oxidizer  mass  ratio  varies  with 
position  across  the  mixing  zone.  The  nitrogen- oxygen  mass  ratio  also 
varies  with  position  in  the  mixing  zone  because  the  nitrogen-oxygen 
elemental  mass  ratio  in  the  central  stream  is  generally  different  from 
that  for  air.  Rather  than  compute  the  equilibrium  mixing  zone  tem¬ 
peratures  and  gas  properties  with  the  usual  cumbersome  equilibrium 
chemistry  calculation  procedure,  a  technique  has  been  developed,  in 
which  the  tabulated  equilibrium  properties  for  fuel -air  mixtures 
(e.  g. ,  Refs.  12  and  13)  can  be  used  to  calculate  the  equilibrium  proper¬ 
ties  of  the  mixing  zone  gases  which  have  variable  nitrogen- oxygen 
mass  ratio. 

NOMENCLATURE 

The  nomenclature  used  in  this  appendix  is  as  follows; 

cp  Specific  heat 

gn  Conversion  constant 

h  Static  enthalpy  of  reactants  (per  unit  mass) 

J  Joule’s  constant 

MW  Mixture  molecular  weight 

m;  Mass  fraction  of  the  ith  species 

ni  Mole  fraction  of  the  ith  species  in  reactants 

Nj  Mole  fraction  of  the  i*  species  in  products 

P  Static  pressure  of  mixture 

R  Mixture  gas  constant 

T  Static  temperature 

T0  Stagnation  temperature 

u  Local  flow  velocity 

Unbarred  symbols  refer  to  the  mixing  zone  system  having  a  given 
fuel -oxygen  mass  ratio  and  having  a  given  nitrogen -oxygen  mass  ratio. 
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Barred  symbols  refer  to  the  equivalent  fuel-air  mixture  having  the 
nitrogen- oxygen  mass  ratio  of  air. 


PROCEDURE 


The  General  Electric  tables  for  combustion  gases  (Refs.  12  and  13 
for  hydrogen  and  hydrocarbon  fuels,  respectively)  contain  the  following 
parameters  tabulated  against  T,  for  various  values  of  static  pressure, 
p,  and  fuel-air  equivalence  ratio: 

1.  1  -  Enthalpy  of  the  reactants 

2.  MW  -  Molecular  weight  of  products 

3.  Nh2  ,  No2  -  Mole  fractions  of  products 

Consider  one  mole  of  products  of  the  equivalent  fuel-air  mixture  at 
a  given  equivalence  ratio.  At  any  temperature,  T,  and  static  pressure, 
p,  the  equivalent  fuel-air  system  is  completely  defined.  If  the  nitrogen 
is  assumed  to  be  inert,  then  the  only  effect  of  the  nitrogen  on  the  equilib¬ 
rium  chemistry  of  the  reacting  species  is  in  the  relationship  of  the  species 
partial  pressures  to  the  mixture  pressure. 


In  the  fuel-air  system 

(  nN2 /n0j  )  =  3‘76 

In  the  variable  nitrogen  system,  the  nitrogen-oxygen  ratio  is  deter¬ 
mined  from  the  mixing  profiles  as  a  mass  ratio,  mN2/mo2 ,  and 

nN2/no2  =  1-142  mN2/m02 

By  considering  one  mole  of  fuel -air  products,  a  certain  fraction  of  the 
nitrogen  must  be  removed  to  convert  to  the  variable  nitrogen  system. 


A  NNj  =  Nn2  -  NNj  =  NN 


(i  _  nw»/nM 

J  ^  3.76  ) 


(I-D 


where  Nn2  is  the  nitrogen  mole  fraction  in  the  fuel-air  products.  The 
resulting  number  of  moles  is  (1  -  A  Nn2),  and  the  mass  removed  is 
28.  016  ANn2  •  The  molecular  weight  in  the  mixture  is  then 

MW  -  28.016  AnNj 


MW  =  . 

I  -  ANNj 

The  resulting  mixture  pressure  is 

P  =  p  (1  -  ANNj> 

The  enthalpy  of  the  reactants  in  the  variable  nitrogen  system  is 

h  = 


MW  h  -  (28.016  ANXl)  hNj 


MW  -  28.016  ANn 


(1-2) 


(1-3) 


(1-4) 
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The  mole  fractions  of  the  products  (except  nitrogen)  are 
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The  mole  fraction  of  nitrogen  in  the  products  is 

,  (1-6) 

The  calculation  procedure  used  to  determine  the  mixture  temperature, 
gas  properties,  and  composition  of  an  equilibrium  system  with  mixture 
pressure,  p  ,  and  reactant  enthalpy,  h  ,  is  as  follows: 


1.  Assume  an  equivalent  system  pressure,  p  . 

2.  Guess  the  static  temperature  of  the  products,  T  . 

3.  Calculate  AK\a  ,  MV  ,  p  ,  h  ,  and  product  mole 
fractions  from  Eqs.  (1-1)  through  (1-6). 

4.  Repeat  steps  (2)  and  (3)  until  the  calculated  value  of 
h. equals  the  initial  value,  thus  defining  T. 

5.  Repeat  steps  (1)  to  (4)  until  the  calculated  value  of  p 
equals  the  initial  value.  The  results  vary  slowly  with 
p  ;  therefore,  a  satisfactory  procedure  is  to  choose 
values  of  p  so  that  the  calculated  values_of  p  bracket 
the  initial  value.  The  results  for  h  ,  MW  ,  and  mole 
fractions  are  then  obtained  by  interpolation. 


MIXTURE  PROPERTIES 

The  mixture  specific  heat  is  evaluated  on  a  mass  basis  as 

c p  —  i  m  j  c 

The  mixture  gas  constant,  on  a  mass  basis,  is 

R  =  BTU  /  lbm  °R 

The  mixture,  y  ,  is  given  by 


y  = 


-  R 


For  a  given  local  flow  velocity,  u  ,  the  local  stagnation  temperature, 
T„  ,  is  evaluated  from  energy  equation: 

T0  =  T  +  - h! - 

^  cp  Sc  J 


19 


AEDC-TR-65-4 


REFERENCES 


1.  Pool,  H.  L.  and  Charyk,  J.  V.  "Theoretical  and  Experimental 

Investigations  of  the  Mixing  of  a  Supersonic  Stream  with  an 
Induced  Secondary  Stream  as  Applied  to  Ducted  Propulsive 
Devices.'1  Project  Squid  Technical  Report  No.  25,  Princeton 
University,  September  1950. 

2.  Perini,  R.  L.  ,  Walker,  R.  E. ,  and  Dugger,  G.  L.  "Preliminary 

Study  of  Air  Augmentation  of  Rocket  Thrust.  "  Applied  Physics 
Laboratory  Report  TG  545,  January  1964. 

3.  Libby,  P.  A.  "A  Theoretical  Analysis  of  the  Turbulent  Mixing  of 

Reactive  Gases  with  Application  to  the  Supersonic  Combustion 
of  Hydrogen.  "  General  Applied  Science  Laboratories  Technical 
Report  No.  242,  June  1961, 

4.  Vasiliu,  J.  "Determination  of  Temperature,  Velocity,  and  Concentra¬ 

tion  Profiles  in  the  Mixing  Layer  Between  a  Rocket  Exhaust  Jet 
and  the  Surrounding  Supersonic  Air  Stream.  "  Convair-Astronautics 
Report  ERR -AN-005,  March  1960. 

5.  Mikhail,  S.  "Mixing  of  Coaxial  Streams  Inside  Closed  Conduit.  " 

Journal  of  Mechanical  Engineering  Science,  Vol.  2,  No.  1, 

March  1960,  pp  59-68. 

6.  Yakovlevskiy,  O.  V.  "The  Mixing  of  Jets  in  a  Channel  with  Variable 

Cross  Section.  "  Izvestiya  Akademii  NaukSSSR,  OTN,  Mekhanika 
i  Mashinostroyeniye,  No.  1,  1962,  pp  66-72  (English  translation 
FTD-TT-62  - 1571). 

7.  Abramovich,  G.  N.  "Mixing  Turbulent  Jets  of  Different  Density,  " 

Izvestiya  Akademii  Nauk  SSSR,  Otdeleniye  Tekhnicheskikh  Nauk, 
Mekhanika  i  Mashinostroyeniye,  No.  3,  1961,  pp  55-57. 

(English  translation  FTD-TT-63 -222). 

8.  Abramovich,  G.  N.  "Turbulent  Jets  Theory. "  Gosudarstvennoye 

Izdatel1  stvo  Fiziko  -  Matematicheskoy  Leteratury,  Moscow, 

1960.  (English  translation  -  WP-AFB  Tech.  Documents 
Liaison  Office  MCL-  1256/ 1+2). 

9.  Phares,  W.  J.  and  Loper,  F.  C.  "A  Technique  for  Solving  Integro- 

Differential  Equations  witn  Application  to  Turbulent  Mixing.  " 
AEDC-TDR-64 -209,  November  1964. 

10.  Smith,  R.  E. ,  Jr.,  and  Mat  z,  R.  J.  "Verification  of  a  Theoretical 
Method  of  Determining  Discharge  Coefficients  for  Venturis 
Operating  at  Critical  Flow  Conditions."  AEDC-TR-61-8 
(AD263714),  September  1961. 


20 


AEDC-TR-65-4 


11.  Rhodes,  R.  P.  ,  Rubins,  P.  M.  ,  and  Chriss,  D.  E.  "The  Effect 

of  Heat  Release  on  the  Flow  Parameters  in  Shock -Induced 
Combustion.  "  AEDC-TDR-62-78  (AD275366),  May  1962. 

12.  Browne,  w.  G.  and  Warlick,  D.  L.  "Properties  of  Combustion 

Gases  -  System  H2~Air.  "  General  Electric  Co.,  Cincinnati, 
Ohio,  R62FPD-366,  November  1962. 

13.  H.  A.  Fremont,  et  al.  "Properties  of  Combustion  Gases  -  System 

CnH2n  “  Air .  11  General  Electric  Co. ,  Cincinnati,  Ohio, 

1955,  (2  Volumes). 


to 


Fig.  2  Schematic  of  Duct  Flow 


p*S9-y  ±OQ3V 


ture, 


AEDC-TR-65-4 


Mass  Fraction  of  Central  Stream  Elements,  C - 

104434 

Fig.  4  Mixing  Zon*  Total  Temperature  Function 


Fig.  5  Mixing  2 


Frozen 

Chemistry 


-  F^  Rocket 
O/F  =  3. 2 


1  0.6  0.8 

Df  Central  Stream  Elements,  C 


1 


104435 


Specific  Heat  Function 


w 

o 


Fig.  8  Schematic  of  Experimental  Apparatus 


*-S9-H10a3V 


a.  Rake  Geometry 


All  Dimensions  in  Inches 
-Rubber  "0"  Ring 


Water  Manifold 


VVVWT 


Pitol 


b.  Probe  Tip  Detail 
Fig.  9  Survey  Rake 


Sapphire  Watch  Bearing 


oo 

CO 


Fig.  10  Apparatus  far  Incompressible  Experiments 
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Fig-  12  Centerline  Velocity  for  Incompressible  Mixi 
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a.  Schematic  of  Initial  Rocket  Pluming  Effect 


Treated  as  Initial 
Section  for  Mixing 


b.  Approximate  Treatment  of  Rocket  Pluming 
Fig.  13  Initial  Pluming  of  Racket  Jet 
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TABLE  I 

EXPERIMENTAL  PARAMETERS  FOR  MIXING  AND  BURNING  APPARATUS 


Nominal  Rocket  Parameters 


Oxidizer 

Gaseous  Oxygen 

Fuel 

Gaseous  Hydrogen 

Chamber  Pressure  (p0j) 

315  psia 

Calculated  Vacuum  Thrust 
(Fj) 

391  lbf 

Characteristic  Velocity 

(c*) 

8200  ft/sec 

Mixture  Ratio,  by  mass 
(0/F) 

3.  16  (±3  percent) 

Total  Propellant  Mass 

Flow  (wj) 

0.  97  lbm/sec 

Throat  Diameter 

1 .  00  in. 

Nozzle  Exit  Radius  ( rNE ) 

1.  19  in. 

Nozzle  Configuration 

15-deg;  half-angle, 
conical 

Mixing  Duct 

Inlet  Diameter  (Dt) 

6  in. 

Length  (L) 

15.  5  in. 

Divergence  half-angle 

3  deg 

Nominal  Air  Temperature  (Toa) 

530°R 
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